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Abstract 

A phenomenological model of parametric surface waves (Faraday waves) is introduced in 
the limit of small viscous dissipation that accounts for the coupling between surface motion 
and slowly varying streaming and large scale flows (mean flow). The primary bifurcation of the 
model is to a set of standing waves (stripes, given the functional form of the model nonlinearities 
chosen here) . Our results for the secondary instabilities of the primary wave show that the mean 
flow leads to a weak destabilization of the base state against Eckhaus and Transverse Amplitude 
Modulation instabilities, and introduces a new longitudinal oscillatory instability which is absent 
without the coupling. We compare our results with recent one dimensional amplitude equations 
for this system systematically derived from the governing hydrodynamic equations. 
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I. INTRODUCTION 



The purpose of this paper is to couple a phenomenological order parameter model 
of parametric surface waves in the limit of weak viscous dissipation to slowly varying 
flows (mean flows). To date, most theories of parametric surface waves near onset have 
neglected such flows despite the observation that their effect is of the same order as other 
cubic nonlinear conservative terms retained. The coupling to the phenomenological model 
presented here allows us to discuss the simplest consequences that these flows have in a 
laterally unbounded geometry, namely shifts in thresholds of secondary instabilities of 
the base pattern of standing waves, and the appearance of a new longitudinal oscillatory 
instability. 

When a layer of an incompressible fluid is vibrated periodically along the direction 
normal to the free surface at rest,it can exhibit parametrically driven surface waves, also 
known as Faraday waves Q, Just above the primary instability of the planar free 

surface, a set of standing surface waves emerge leading to a stationary pattern with a 
symmetry that depends on the physical parameters of the fluid and the frequency of the 
forcing ji], Is , E , 111 ■ Intricate phenomena appear in the limit of weak viscous dissipation in 
which nonlinear wave interactions responsible for wave saturation and pattern selection 
are dominated by triad resonant interactions ^| Q, \v\ . Whereas the first bifurcation 
away from planarity is to a set of standing waves in which mean flow effects are absent, 
mean flows are expected to be important in determining the stability of the primary 
waves, and more generally in weakly damped systems. In this latter case, standing wave 
amplitude equations can be expected to be valid only very close to onset. 

Current weakly nonlinear theory is restricted to the small region above threshold in 
which standing waves are stable, a state in which mean flows identically vanish. However, 
the contribution from mean flows to the equations governing the slow evolution of the 
surface waves can be of the same order as the standard cubic nonlinear and conservative 
terms which are always retained. Thus weakly nonlinear corrections to surface waves 
and mean flows must be considered simultaneously, which has not been done in three 
dimensions in the limit considered below (see for the analysis of this limit in two 

nnn 

dimensions and [13, 14L llq for the analysis of related limits). The effects considered here 
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are not unlike other kno wn p henomenology that includes the streaming flow produced by 
fixed surface waves, [lfl llT], llSl and references therein, and the evolution of surface 
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waves in the presence of a fixed vortical flow 

A consistent introduction of mean flow effects into the amplitude equations for Faraday 
waves requires explicit consideration of special limits that involve the physical dimensions 
of the container. We specifically focus here on the case of a fluid depth that is logarith- 
mically large compared to the wavelength, and derive a set of evolution equations for the 
surface waves and the associated mean flow in the double limit of small viscosity and 
large aspect ratio (the ratio between the lateral size of the container and the wavelength). 
We find two separate contributions to mean flow, namely an inviscid contribution arising 
from the slowly varying motion of the free surface, similar to the one appearing in classical 
Davey-Stewartson models (3], and a viscous one resulting from a slowly varying shear 
stress produced by the oscillatory boundary layer attached to the free surface. This latter 
contribution describes vorticity transport (by viscous diffusion or convection) from the 



boundary layer into the bulk 
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An important simplification in our analysis is that the cubic nonlinear terms of the 
phenomenological model are chosen so as to lead to a stripe pattern above onset instead of 
a square pattern as experimentally observed in the limit of weak viscous dissipation. While 
it is a simple matter to modify the functional form of the cubic term to produce square 
patterns, we have chosen to first clarify the effect of mean flows on slow modulations 
of a stripe pattern. There is no satisfactory theory at present that can account for 
the interaction between slow spatial modulation of the waves and mean flows in three 
dimensions, and the case of stripes is considerably simpler than other symmetries involving 
a larger number of plane wave components at onset. 

II. FORMULATION 

We consider a fluid layer of unperturbed depth d* supported by a horizontal plate that 
is vibrating vertically with an amplitude a* and a frequency 2u*, where the superscript 
* denotes dimensional quantities. In order to facilitate comparison with related results 
in we use for adimensionalization the characteristic time uo* 1 and length k* 1 , 
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where the wavenumber k* is related to uj* by the inviscid dispersion relation 

uj* 2 = g *k* + a*k* 3 /p*, (1) 

in terms of the gravitational acceleration g*, the surface tension a* and the density p*, all 
assumed constant. Here we are assuming that the wavelength k*~ l is (at least, somewhat) 
small compared with the depth of the container. The resulting dimensionless continuity 
and Navier-Stokes equations in a reference frame attached to the vibrating container, with 
the z = plane at the unperturbed free surface, are 

V • u + d z w = 0, (2) 
d t u - w(Vw - d z u) - u^V • u L = - Vp + 7( V 2 u + d 2 zz u) /2, (3) 
d t w + u • (V«) - d z u) = -d z p + 7( V 2 w + d 2 zz w) /2, (4) 

in — d < z < h(x,y,t), with boundary conditions resulting from no slip at the supporting 
plate, 

u = 0, w = at z — —d, (5) 

and kinematic compatibility and equilibrium of tangential and normal stresses at the free 
surface, 

d t h + u • Vh = w, (6) 
d z u + Vw- (Vtt + V« T ) • Vh + [2d z w - (d z u + Vw) • Vh]Vh = 0, (7) 
p - (\u\ 2 + w 2 )/2 - [4asin2t + 1 — T]h + TV • [V/i/(l + |V/i| 2 ) 1/2 ] (8) 
= -f[d z w - {d z u + Vw) • Vh + (V/i • (Vtt + Vu T )/2) • V/i]/(l + | Vh\ 2 ) 

dX z = h. Here 

u = (u,v,0) (9) 
and w are the horizontal and vertical velocity components, 

V=(d X ,dy,0) (10) 

denotes the horizontal gradient, the superscript _L over a horizontal vector denotes the 
result of rotating the vector 90° counterclockwise, namely 

«- L = (-u,«,0), (11) 



and the superscript T over a tensor denotes the transpose; p (=pressure+(|'u| 2 + w 2 )/2 + 
[1 — T + 4a sin(2t)]z) is a conveniently modified pressure, and h is the (vertical) free surface 
deflection. For simplicity we do not consider lateral walls, but impose periodic boundary 
conditions in two horizontal directions, namely 

(u,w,p)(x + L u y,z,t) = (u,w,p)(x,y + L 2 ,z,t) = (u,w,p)(x,y, z,t), 

(12) 

f(x + Lx, y, t) = f(x, y + L 2 , t) = f(x, y, t). 
And for convenience we also consider the vertically integrated continuity equation 



obtained upon integration of (j2J) in — d < z < h and substitution of (|B|) . 
A. Multiple scale analysis: oscillatory and mean flows 

We consider next a specific range of parameters in which it is possible to simplify the 
problem by separating fluid motion into a "fast" oscillatory component, and a "slow" 
mean flow. In particular, we consider the system of surface waves near onset, and in the 
limits of a very large lateral surface and weak viscous dissipation. The problem depends 
on the following dimensionless parameters: the dimensionless viscosity 
(with v* the kinematic viscosity), the gravity-capillary contribution T = a*k* 3 /(p*u* 2 ), 
the forcing amplitude a = a*k*, the container depth d = d*k* and the aspect ratios L\ and 
L 2 ; note that, according to (JTJ, < V < 1 and the extreme cases T = and 1 correspond 
to the purely gravitational and purely capillary limits respectively. The approximation 
below requires that (a) the aspect ratio of the container be large, (b) the surface waves 
be weakly damped and (c) exhibit a small wavelength compared to the container's depth 
and (d) a small steepness, which in turn require that 



where L < min-fLx, L 2 }. The large spatial scale set by the (large) aspect ratio introduces 
a slow horizontal scale over which both spatial and temporal modulations are expected 
to occur. As suggested by the 2-D case |l2|, this scale is expected to be determined (in 




(13) 



L > 1, d> 1, 7 < 1 



V/i| < 1, a < 1, 



(14) 
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the equations for the oscillatory flow associated with the surface waves) by the balance 
between cubic nonlinearity and either (i) transport with the group velocity or (ii) disper- 
sion. And d must be not too large, see below. For the sake of clarity we assume that 
d is logarithmically large compared to the remaining small parameters (namely 7, a and 
and we shall treat d as a 0(1) parameter. In fact, for simplicity we consider the 
distinguished limit 

7 2 < e~ d < L' 1 ~ 7 ~ a ~ e 2 , (15) 

where e is a measure of the surface wave amplitude, see (j!8)) below. The estimates 
7 ~ a ~ e 2 result from imposing that linear damping, cubic nonlinearity and parametric 
forcing be of the same order. Therefore we are implicitly assuming that the coefficient 
of the cubic term is 0(1), which excludes triad resonances [llj]. If these are present, the 
cubic coefficient becomes 0(7 _1 ) and a different scaling applies. In order to concentrate 
on the effects of mean flows, we exclude triad resonances in what follows. 

Under these assumptions we shall (implicitly) use a multi-scale analysis in both (hor- 
izontal) space and time. But in order to make a not too technical presentation and to 
avoid obscuring the main ideas with a too involved notation, we shall use only one time 
variable and one space variable in each horizontal direction. The basic (fast) scales involve 
0(1) increments of x, y or t. A magnitude ip that exhibits these fast scales at leading 
order is such that 

\dip/dx\ ~ \dip/dy\ ~ |V| or \dip/dt\ ~ (16) 

and it will be said to depend strongly on the associated variable x, y or t. If instead the 
magnitude only changes over the slower scale at leading order, namely if 

\dip/dx\ < M, \dif>/dy\ < H or \di/>/dt\ < (17) 

then the magnitude will be said to depend weakly on the associated variable, x, y or t. To 
proceed, we decompose the flow variables and the free surface deflection into oscillatory 
and time-averaged parts, associated with the surface waves and the mean flow (denoted 
hereinafter with the superscripts o and m), respectively, as 

(u,w,p,h) = e(u ,w ,p ,h°)+e 2 (u m ,w m ,p m ,h m ), (18) 
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where (i) the (oscillatory) flow variables associated with the surface waves are required to 
be such that 

(u°) ts = 0, (w°) ts = (p°) ts = (h°) ts = 0, (19) 
with (-) ts standing here and hereafter for the time average in the basic oscillating period 

rt+2n 

(^) ts = (2n)- 1 if;dt; (20) 



and (ii) the variables associated with the mean flow are required to depend weakly on 
time; more precisely we assume that 

\d t u m \ ~ e 2 \u m \, \d t w m \ ~e 2 \w m \, \d t p m \ ~ e 2 \p m \, \d t h m \ ~ e 2 \h m \, (21) 

where we are anticipating the timescale for the slow evolution of the mean flow, t ~ e~ 2 . 
Also we anticipate that the re-scaled flow variables u°, . . . , h°, u m , . . . , h m are at most of 
order unity; see below. The mean flow is described in terms of the time-averaged velocity 
e 2 u m , which is the Eulerian velocity and does not coincide in general with the velocity 
associated with the time average of the trajectories of fluid elements. The latter is the 
Lagrangian mean velocity, or mass transport velocity (denoted here as e 2 u mt ) which is 
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the appropriate velocity to analyze mean trajectories of passive scalars (see e.g., 
in connection with chaotic advection j^j]). The difference between the two is the Stokes 
drift so that, in the notation of this paper, its horizontal and vertical components, scaled 
with e 2 , are given by [3] 



u mt - u m = u Sd = ((J V<f)-V)V(j)+(J d z <f))d z V(j)} ts 
((J V<p)-V(d z <P) + (J dA)d zz $) 



, (22) 

w mt _ w m — w Sd — III T-7A\ . T-TftS A\ I ( I a A\ts 

in a first approximation, where we are anticipating Eq. (|2*K|) below, and the operator J* 
is defined as 

$={( ^dt) to , (23) 

J t 

with (-)* standing hereinafter for the time-oscillatory part, defined as 

(^) to = ^) ts . (24) 
By definition, Eq. (|2~H|) is independent of to- 



B. Oscillatory flow 



We begin by deriving the equations governing the oscillatory flow associated with the 
surface waves which exhibits a thin viscous boundary layer of 0{J^f) thickness attached 
to the free surface. In the bulk region outside this boundary layer the oscillatory velocity 
components and the pressure are given by 



u° = V<p + e% d z <P)(yw m -d z u m ) + ( / V0)^V • u m± ] + 0(e d 

J r t J (25) 

w° = d z <f)-e 2 ( V<f))-(Vw m -d z u m ) + 0(e 3 ), p=-dt<f>, 

as obtained upon substitution of Eq. (JTHJ) into Eqs. ©-(01), where is the velocity 
potential. Thus the oscillatory flow is potential at leading order, but not at order e 2 
which must be retained in what follows. Substitution of Eq. (|23|) into Eq. (J2J) yields, 
after some algebra, 

r-t r t 

72 i i q2 / i 2/ / a i\T7 /T-7„..m o „.m\ , Jl 



V 2 (f> + d 2 zz (j) + e 2 (J d z( j))V - (Vw m - d z u m ) + e 2 V -[(J V^^-tt" 11 ] 

~ e2 (f V0) ■ d z (Vw rn - d z u m ) =0 for -d<z<0. (26) 

Here we are taking the upper boundary at the unperturbed free surface, which can be 
done because h is small. The boundary conditions at the upper boundary must include 
the effect of the vortical flow in the boundary layer attached to the free surface on the 
oscillating flow in the bulk. To the approximation relevant here this only requires to 
replace the boundary condition by (see, e.g., jll| ) 

d t h + u • Vh = w + W(x, y, t) at z = h, (27) 

where W is given by 

d t W(x,y,t) = 1 V 2 (d z( py, (28) 
this equation can be integrated to obtain 

W = 1 j V 2 («9 2 0). (29) 
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Thus, to the approximation relevant here Eq. (jTHj) can be rewritten as 

d t h + V -(j ud^j= 1 j V 2 (<9,0), (30) 

The boundary conditions for the oscillatory flow at the unperturbed free surface are now 
obtained by a Taylor expansion from Eqs. (JHjK and (j3Uj) . and are found to be 

d t h° - d z <p + e 2 { J' V0) • {Vw m - d z u rn ) + £ (V • {h°V<P)) to 

+ e 2 V ■ [h m V<P + h°u m + [h°) 2 d z V<P/2] = 7 J V 2 {d z <P), (31) 

d t <t> + e(h°d 2 tz <P) to + e 2 [h m d 2 tz <P - h°d zP m + (h°) 2 dl z <P/2\ + e(| V0| 2 + \d z <P\ 2 ) to /2 

+ e 2 [u m • V0 + h°d z (\V<j)\ 2 + \d z <j)\ 2 )/2]+4:a(h sm2t) t0 (32) 

+ (i - v)h° - rv • [v/i°/(i + e 2 \ vh°\ 2 y/ 2 ] + 7 af> = o, 

where we are using Eq.fJo'Hj) below. The boundary condition (jSJ) at the lower plate and 
the periodicity conditions (fT2j) yield 

d z (j) = at z = —d, (33) 
(j)(x + L h y,z,t) = <p(x,y + L 2 ,z,t) = <f>(x,y,z,t), 

(34) 

/i°(x + L 1 ,y,t)^ h°(x, y + L 2 , t) = h°(x, y, t). 

We are consistently neglecting terms of order e 3 in Eqs. (j3*T|) - (jH2^ because of the 
approximations listed in Eq. (|15|) . and taking into account that those terms that are 
either (a) cubic in the oscillatory flow variables or, (b) linear in both a oscillatory variable 
and a slowly varying variable, exhibit zero temporal mean values at leading order. 

Before proceeding any further we note that mean flow does not contribute to the 
averaged (in the time scale t ~ 1) energy equation at leading order, which is consistent 
with the fact that mean flow variables (velocity and free surface deflection) are small 
compared to their counterparts in the oscillatory flow. The averaged energy equation is 
obtained upon multiplication of Eq. (|2T)j) by d t <fi, integration in < x < Li,0 < y < 
L2, —d < z < 0, averaging over a period of oscillation, integration by parts repeatedly 
and substitution of Eqs. (jSU)-© and flM|). We find 

\(d z (\V<P\ 2 + \d z <P\ 2 ) + 8a(d t h )(h°sm2t) t0 ) ts dxdy + 0(-f + a + e 2 ), (35) 



d£_ 

dt 



Jo 
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where £ is the time averaged (kinetic plus potential) energy. The first term in the 
right hand side of Eq. (|H5jl (except for a factor of two) is the classical result, first 
given by Landau and Lifshitz |29| that approximated viscous dissipation by linear damp- 
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111]). Note that mean flow variables 



ing from the bulk potential flow (see also 
(both velocity and free surface deflection) are small compared to their counterparts in 
the oscillatory flow and do not contribute to the energy at leading order. To obtain 
Eq. (|H5jl we have taken into account that u m and w m are independent of t at lead- 
ing order, and that if ip and (p are t-periodic, of period 2ir (as the variables associ- 
ated with the oscillatory flow are to first approximation), then to leading order we have 
((p f ip) ts = -(tpf y) ts and (ip f tp) ts = 0. 

C. Mean flow 

In order to obtain the equations and boundary conditions governing the slowly varying 
flow we must take into account the oscillatory boundary layer attached to the free surface, 
which provides (at the edge of this layer) a slowly varying shear stress that must be 
imposed as a boundary condition for the mean flow in the bulk. This forcing mechanism 



was first uncovered by Longuet-Higgins 



23| . who obtained an explicit expression for 



the forcing shear stress produced by general boundary layers in 2-D. The counterpart 
of this expression in 3-D (for a free surface of general shape) has been only obtained 
quite recently (^J, although quasi-planar free surfaces (as the ones considered here) were 
already considered in a not too well known work [3J|. With the notation of this paper, 
the general formulae derived in j^O] yield 

<9 2 u m + Vw m = 2(V (V • (/i°V0)) + (V/i° • V)V<f)+{V 2 <f))Vh ) ts at 2 = 0, (36) 

where only the leading order contribution as 7 — > and e — > is retained. The boundary 
layer attached to the free surface has no effect on the other two boundary conditions at 
the unperturbed free surface, which are directly obtained from Eqs. (JEJ) and |T3|) to be 



p m - (1 - T)h m + rV 2 F = (h dl<P + (| V0| 2 + |<9 2 0| 2 )/2)* s (37) 

and d t h m + V • ( / u m dz) = - V • ((/i°V0)' s ) at z = 0, (38) 
J-d 
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where we are only taking into account the leading order terms. And from Eq. (|SJ) we have 
no slip at the lower plate 



u m = 0, w m = at z = -d, (39) 

at leading order. We are neglecting the effect of the oscillatory boundary layer attached 
to the lower plate because its effect is quite small (the horizontal component of the mean 
flow velocity near the lower plate is proportional to the square of the vertical jump of 
the horizontal component of the oscillatory velocity accross the lower boundary layer 2^ |. 
which is 0(e~ 2d )e 2 <C e 2 (Eq. p4jl ): this in turn is small compared to the streaming flow 
velocity in the bulk, which is 0(e 2 )). These boundary conditions show that mean flow is 
forced by surface waves in two ways. Those terms appearing in the right hand sides of 
Eqs. (p?Tj) and (jHHj) provide an inviscid forcing mechanism that by itself would provide an 
inviscid mean flow, like that appearing in the Davey-Stewartson model 2^]. The right 
hand side of Eq. (}36|) instead produces a forcing shear stress that drives a viscous mean 
flow, which is absent in the usual inviscid and nearly inviscid theories of Faraday waves. 
Note that his forcing stress is generically non zero and independent of viscosity at leading 
order, fact that is well known but somewhat surprising because this effect is due to the 
oscillatory boundary layer, and is absent in the strictly inviscid case. We remind the 
reader that the limit of vanishing viscosity is a singular limit which does commute with 
the limit e — » 0. We could decompose the mean flow into its inviscid and viscous parts, 

nn 

as is done in |13l . Il4j , but for convenience this is not done here. 

Finally, we substitute Eqs. (fT£|) and (f23|) into Eqs. and 0, and take the time 
average defined in Eq. (J2*U|) in the resulting equations. Proceeding as we did to obtain 
Eq. we find, after some algebra 

V • u m + d z w m = 0, (40) 
d t u m - e 2 [{w Sd + w m ){Vw m - d z u m ) + (u Sd + u m ) L V ■ u m± ] 

= -Vq m + 7 (V 2 <u m + d 2 zz u m )/2, (41) 

d t w m + e 2 {u Sd + u m ) ■ {Vw m - d z u m ) = -d z q m + 7 (V 2 w m + d 2 zz w m )/2, (42) 

where u Sd and w are the horizontal and vertical components of the Stokes drift given 



11 



by Eq. (J22J), and the modified pressure q m is defined as 

q m =p m + e 2 (d z (l)(J V0) • (Vw m - d z u m )) ts . (43) 
Finally the periodicity condition (|T2*|) yields 

(u m , w m , q m )(x + L u y, z, t) = (u m , w m , q m )(x, y + L 2 , z, t) = (u m , w m , q m )(x, y, z, t), 
h m (x + L 1 ,y,t) = h m (x, y + L 2 , t) = h m (x, y, t). 

(44) 

In order to estimate the magnitude of the various terms that depend on the surface 
wave variables and force the mean flow, we assume that at leading order the velocity 
potential and the free surface deflection h° can be written as a superposition of plane 

waves 

N N 

= ie 2 A n e i{t+kn ' x) + c.c. + . . . , h° = ^ A n e l{t+kn ' x) + c.c. + . . . , (45) 

n=—N n=-N 

where the complex amplitudes, A_jy,. . . ,An are allowed to depend only on slow space and 
time variables. The wavevectors k_^,. . . ,kjy correspond to only iV directions because they 
are related in pairs as 

k_ n = —k n and |fc n | = 1 forn = 1, . . . , N. (46) 

Thus for each n = 1, . . . , N the n-th and the (— n)-th waves counter-propagate along the 
same direction. Note that each pair of counter-propagating waves builds a standing wave 
in the short timescale if and only if \A n \ = |A_ n |, and the whole surface wave pattern 
is seen to be standing in the short time scale if and only if the following, more stringent 
condition holds 

A m A n = A_ m A_ n for all m, n = 1, . . . , JV. (47) 

Such surface waves will be called quasi- standing below. 

By using Eq. (J45|) . the forcing terms in the boundary conditions Eqs. (l36|) - (}38|) are 
written as 

(2h°dl<P + | V0| 2 + \d z <P\ 2 ) ts = (| V0| 2 - \d z <P\ 2 ) ts = J2(k m ■ k n - l)A m A n e^~ k -> x 

m,n 
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+ KAiK • VA m - A m k m • VA,)^-*" 1 " + c.c. + 0(L~ 2 ), (48) 

m,n 

V • (/i°V0) is = (V • h°V(f>) ts = -2 J> m • V(\A m \ 2 ) + 0(L- 2 ), (49) 

m 

(V (V • (h°V(j>)) + (Vh° • V) + (V 2 0) V/i°) ts = 

^(1 + k m • fc„,)fc m A m A n e i(fc "- fc '" ) - :c + c.c. + 0(L _1 ), (50) 

at 2 = 0; and the Stokes drift velocity components in Eq. fl22(l are written as 

u s d = w -i e 2, ^- (fcm . fcn + i) km A m A n e^- k -> x + c.c. + OiL- 1 ), (51) 

w 5d = 2uT 1 e 2 * ^ fc m • V(|A m | 2 ) + 0(L- 2 ), (52) 

m 

in —d < z < 0. 

Several remarks can be made about these boundary conditions. First, the forcing term 
given in Eq. f)48|) depends only weakly on the horizontal variables x and y, namely 

+ I V0| 2 + \d z cl>\ % ) ts ) hs ~ IT 1 , (53) 
where is the horizontal average in the short spatial scales 

(i/j) hs ^(J dxdy) J ipdxdy. (54) 

Here B is a ball of radius large compared to 1 but small compared to L. Or, in terms of 
a horizontal Fourier transform with associated wavenumbers k mn ^ if (m, n) 7^ (0,0), 
with fc o = 0, 

(tf>) ha = ^oo if V> = X^e 1 *" 1 """, (55) 

m,n 

where ip mn is allowed to depend weakly on x, y and t (and strongly on z). 

Second, the forcing terms in Eqs. (J4*9~j) and (J5U|) . and the Stokes drift vanish at leading 
order if the surface wave pattern is quasi-standing 

(V (V • (h°V(j))) + (V/i° • V)V0 + (y 2 <p)Vh°) ts ~ w Sd ~ L" 1 , 

and V • (/i°V0) is ~ w 5d ~ L" 2 if Eq. flSJ) holds. (56) 
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Third, and according to Eqs. (f47Ij) and (JH2*|) . we have 

| [(V • (^°V0)) is ] 2=o | = OiL- 1 ), \w Sd \ = OiL- 1 ). (57) 

By using the continuity equation (|40|). the boundary condition (jHH|) can be rewritten as 
9 t /i m + w m + (/i°V0)* s ) = 0, which invoking Eqs. (|2Tj) and (JH7|) yields, at leading order, 

w m ~ L" 1 at 2 = 0. (58) 
i. Short and long wave decomposition of the mean flow 

We next decompose the mean flow variables into a short wave component (or oscillatory 
in the horizontal direction) and a long wave component (or slowly varying in the horizontal 
direction) 

(u m , w m , q m , h m ) = (u mo , w mo , e 2 q mo , h mo ) + (u ms , e 2 w ms , q ms , h ms ), (59) 

where the short wave component is such that 

(u mo ) hs = 0, (w mo ) hs = (q mo ) hs = (h mo ) hs = 0. (60) 

The long wave components depend weakly on the horizontal variables 

|Vu ms | ~ £ 2 \u ms \, \Vw ms \ ~e 2 \w ms \, \ Vp ms \ ~ e 2 \p m % \ Vh ms \ ~ e 2 \h ms \, (61) 

where we are assuming that the slow spatial scale for horizontal gradients of the long 
wave mean flow is the same as that of the envelope of surface waves, namely of the order 
of L ~ e- 2 (see flEJ). Also, in Eq. (H3J) u mo ~ w mo ~ q mo ~ h mo ~ u ms ~ w ms ~ g ms ~ 
fi ms ~ 1 and thus we are anticipating the order of magnitude of all variables associated 
with both mean flow components. 

The equations governing the short wave component of the mean flow are obtained by 
substituting Eq. (J5T?|) into (|3fi|l - (|4"2|l . The short wave deflection h mo is given by 

_ (i _ T)h™ + TV 2 h m ° = [((h°dU + (I V0| 2 + \d z( j ) \ 2 )/2) ts ) ho ] z=0 , (62) 

with periodic boundary conditions h°(x + Li,y,t) = h°(x,y + L 2 ,t) = h°(x,y,t), and 
where where {■) ho denotes the short wave component 

W h ° = V - W hs - (63) 
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Note that, according to Eq. (JHHj) . the horizontal mean value of the right hand side of Eq. 
(I62j) vanishes at leading order, as required by volume conservation. 
Short wave velocity and pressure u mo , w mo and q mo are given by 

V • u mo + d z w mo = 0, (64) 
d t u mo - e 2 (w mo (Vw mo - d z u mo - d z u ms ) + (u Sd + u mo + w ms ) x V • u mo± ) h ° = 

-e 2 Vq mo + 7 (V 2 w mo + d 2 zz u mo )/2, (65) 
d t w mo + e 2 ((u Sd + u mo + u ms ) • {Vw mo - d z u mo - 9 z w ms )) ho = 

- e 2 d z q mo + 7 (V 2 w mo + d 2 zz w mo )/2, (66) 

in — d < z < 0, with boundary conditions 

w mo = 0, w mo = at ^ = -d, (67) 
^w mo = 2(((V/i°- V)V0+(V 2 0)V/i°) ts )' 10 , ^ mo = at 2 = 0, (68) 
(u m °, w mo , q mo ) (x + L u y, z, t) = (u m °, w mo , q mo ) (x, y + L 2 ,z, t) 

= (u mo ,w mo ,q m °)(x,y,z,t), (69) 
h mo {x + L u y,t) = h mo {x, y + L 2 , t) = h mo {x, y, t). 

where we have taken into account Eq. (|57|). Note that this problem is decoupled from 
that giving h mo (Eq. 

In order to determine the equations governing the long wave component of the mean 
flow, we first take into account that, according to the continuity equation and the bound- 
ary condition (JHHJ), the rescaled vertical velocity component is given by 

w ms = _ £ -2 f v . u ms ^ ^ 
J-d 

and, according to Eq. (jHTjl . is of order unity as assumed above. The problem giving 
(u ms , q ms , h ms ) becomes decoupled from w ms as we show now. From the momentum 
equation (|4~2"|) we obtain 

q ms = q ms (x,y,t), (71) 
and then, by invoking Eq. (jBTjl . the momentum equation (|4*T|l yields 

d t u ms -e 2 (w m %Vw mo -d z u mo ) + (u Sd + u mo ) ±J V-u mo± ) hs = -Vq ms + -fd zz u ms /2 (72) 
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in —d < z < 0, with boundary conditions 



u ms = at z = -d, d z u ms = 2(((W- V)V0+ (V»V/i°) ts ) /ls , (73) 

r° 

q ms -(1-T)h ms = 0, d t h ms + V-( u ms dz) = -V • (((h°V <p) ts ) hs ) , (74) 

at z — 0, resulting from Eqs. (|3*B |1 -(|3* ^1 and The periodic boundary conditions (fT2*j) 
lead to 

(u ms , w ms , q ms ) (x + L 1 ,y,z,t) = (u ms , w m \ q ms ) (x, y + L 2 , z, f ) 

= (u ms ,w ms ,q ms )(x,y,z,t), (75) 
/i ms (z + U , y, t ) = h ms (x, y + L 2 , t) = h ms (x, y, t ) . 
All terms appearing in Eqs. (|72 p -(|74j ) are of the same order because 7 ~ e 2 (see Eq. 
IfHijl) and u ms , q ms and h ms satisfy Eqs. (J2IJ) and (jHH). The forcing term in Eq. Ij73h) 
V • ((/i o V0)* s ) /ls , is inviscid, and is the only forcing term that appears in the standard 



Davey-Stewartson model [22j which involves a potential mean flow. This model is only 
valid for stripes patterns as described below in Eqs. (|83 j) - (J55|) . where this term is precisely 
the only forcing term remaining in the right hand side of Eq. (|86|). On the other hand, 
the forcing terms in the right hand sides of Eqs. (}68|) and (J73*|) are necessarily viscous and 
drive a vortical flow. 

Substituting Eqs. (JHHj) into Eqs. (|2l)|) - (J3H|) . we obtain the following equation (after 
some algebra): 



V 2 (j) + d 2 zz (j) + e 2 (J d z (f))V -(Vw mo -d z u mo )+e 2 V -[(J V0) x V • u mo± ] 

-e 2 (J V<p)-d z {Vw mo -d z u mo -d z u ms ) = (76) 
for — d < z < 0, and boundary conditions 

d t h° - d z <f) + e 2 { [ V0) • {Vw mo - d z u mo - d z u ms ) + e{V • {h°V<f))) to 



+ e 2 V • [(h mo + h ms )V(f) + h°(u mo + u ms ) + (h o yd z V<p/2} 

= 1 j t V 2 (d z( p), (77) 

d t <P + e(h°d 2 z <P) to + e 2 [(h mo + h ms )dU - h°d z q ms + (h°) 2 d* z J /2] 

+ e<| V0| 2 + \d z <p\ 2 ) to /2 + e 2 [{u mo + u ms ) ■ + V0| 2 + |^0| 2 )/2] 
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+ 4a(h°sm2t} to + (1 - T)h° - TV • [W/(l + e 2 \ V/i°| 2 ) 1/2 ] + 7 <9 2 > = 0, (78) 
d z (j) = at z = -d. (79) 

This is the central result of this section. In the limit considered (Eq. 1)15)1 ). flow 
variables have been decomposed into oscillatory and slowly varying parts (Eq. ([18))). 
where the oscillatory components are given by Eq. (|25p. with </> and h° given by Eqs. 
(126)) - The slowly varying component has been further decomposed into short wave 
and long wave components in Eq. (|59|) . with the short wave component given by Eq. (|62|) . 
and the long wave given by Eqs. (f?2*|) - (f?l| . Thus the coupled evolution of surface waves 
and mean flow is given by Eqs. (J02J) , flM!)-®, (|72jl - l(75jl . and Ij7^ - (|75)l . This system of 
equations still includes the full 3-D Navier-Stokes equations and (jMj) with a large 
Reynolds number based on the horizontal size (although a 0(1) Reynolds number based 
on the containers depth). 

The numerical solution of this coupled problem remains quite complicated, and further 
simplifications are necessary to make it tractable. We discuss in the following section a 
hierarchy of simplified models that are based on various physical assumptions and in some 
cases on ad hoc approximations. 

III. APPROXIMATIONS TO THE COUPLED MEAN FLOW-SURFACE WAVE 
EQUATIONS 

A. Stripes 

In this particular case the equations derived above simplify substantially. Unfortu- 
nately in the nearly-inviscid limit we are considering stripes are not the selected pattern 
in 3D (3^], except in the purely gravity wave limit of F <C 1. This would require that the 
basic wavelength 2n/k* be large compared to the capillary length £ c = a/<t/ (pgo), which 
is of the order of 3 mm for water. We consider first the strict gravitational limit of T — 0. 
If only one stripe is present at each point then we can take N = 1 in Eqs. ()49)) - (fH2*|) at 
this point, and the short wave parts of u Sd and ((V/i° • V) + (V 2 </>) ^h°) ts vanish in 
Eqs. (!64|) - (p)B*j) . This implies that the short wave part of the mean flow is unforced, and 
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thus admits the solution u mo = 0, w mo = 0, q mo = constant which we assume is globally 
stable. Then for large times such that the short wave component vanishes, the remaining 
variables are given by Eqs. (HJ), (j3Tjl - (j33jl . (O and (f72]l - (j73|l . By also invoking Eq. CTl . 
these latter equations can be rewritten as 

V 2 + <9 2 > + e 2 ( J* V0) • <9^u ms = 0, (80) 
9 iW ms = -V/i ms + jd zz u ms /2 (81) 

in — d < z < 0, with boundary conditions 

<9 Z = w ms = 0, (82) 

at z — —d and 

d t h° - d z <p - e 2 { J V0) • d z u ms + e(V • (h°V<p)) to 

+ e 2 V -[(h mo + h ms )V(j) + h u ms + (h°) 2 d z V(j)/2}= 1 J V 2 (<9 Z 0), (83) 

9,0 + e^WdUr + e 2 [(^° + h ms )dU - d z q ms + {h°) 2 dl z <t>/2\ 

+ e(\ V0| 2 + |^0| 2 ) to /2 + s 2 [u ms ■ V0 + h°d z (\ V0| 2 + |<9 2 0| 2 )/2] 



+ 4a(/i° sin 2t)'° + /i° + 7<9 2 2 2 = 0, (84) 
d z u ms = 2{{{Vh° • V)V0 + (V 2 0) Vh°) ts ) hs , (85) 

d^ ms + V • ( / u ms tfc) = - V • (((h°V<J)} ts } hs ), (86) 



at z = 0. The periodic boundary conditions in the horizontal variables resulting from Eq. 
dI2I) are 

(0, /i°)(x + L 1? j/, z, t) = (0, h°)(x, y + L 2 , z, t) = (0, h°)(x, y, z, t), 

(h°,h mo ,h ms )(x + L u y,z,t) = (h°,h mo ,h ms )(x,y + L 2 ,z,t) = (h°,h mo ,h ms )(x,y, z,t). 

(87) 

These equations are considerably simpler as they only include the heat equation, Eq. 
(jHH), instead of the full continuity and Navier-Stokes equations, but yet allow significant 
variation in stripe orientation. 
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B. Linear approximation for the mean flow 

The mean flow equations and boundary conditions Eqs. (|M j) -([S8 )l and (|72 ]) -(f71 j) are 
linear in the mean flow variables except for convective terms. If these are neglected then 
the mean flow equations are rewritten as Eq. (jo^j) and 



u mo± ] 



V 2 <f) + d 2 zz <p + e 2 {J d z <f))V- (Vw mo - d z u mo ) + e 2 V ■ [{J V^V 

-e 2 (j V<P) ■ d z {Vw mo - d z u mo - d z ii ms ) = 0, (88) 

V • u mo + d z w mo = 0, (89) 

d t u mo = -e 2 Vg mo + 7 (V 2 w mo + d 2 zz u mo )/2, (90) 

d t w mo = -e 2 d z q mo + 7 (V 2 w mo + d 2 zz w mo )/2, (91) 

d t u ms = -Vq ms + -fd zz u ms /2 in - d < z < 0, (92) 

in — d < z < 0, with boundary conditions 

d z <p = 0, u mo = u rns = 0, w mo = w ms = 0, (93) 



at z = —d, and 



d t h° -d z <f) + s 2 (J V0) • (Vw™ - d z u mo - d z u ms ) + e(V • (h°V<P)) to 

+ e 2 V ■ \{h mo + h ms ) V0 + h°(u mo + u ms ) + (h°) 2 d z V0/2] 

= 7 yV 2 (^0), (94) 

d t <p + e{h°d 2 z <P) to + e 2 [(h m ° + h ms )dU - h°d z q ms + (h°) 2 dl z <p/2] 
+ e(\V<f)\ 2 + \d z <f)\ 2 ) to /2 + e 2 [{u mo + u ms ) • 
+ h°d z {\V<f)\ 2 + \d z <f)\ 2 )/2] + Aa(h° sin2t) to 

+ (1 - T)h° - TV • [W/(l + e 2 \ Vh°\ 2 ) 1/2 ] + 7 <9 2 z = 0, (95) 



w 



■ mo - 0, d z u mo = 2(((Vh° • V)V0 + (V 2 <f))Vh } ts } h0 } (96) 



d z u ms = 2(((V/i° • V)V0 + (V 2 0) Vh°) ts ) hs , (97) 

f-0 



(98) 



/■O 

q ms ^ms = Qj ^m. + V • ( / « mfl tfc) = - V • ( ( (h° V<j)} tS } ks ) , 

- (1 - T)h mo + rV 2 /i™ = ((/i°<9 2 z + (| V0| 2 + l^l 2 )^) 45 )' 10 , (99) 
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at z = 0; and periodic boundary conditions in the horizontal direction as in Eqs. ([finj) and 

This linear approximation exactly provides the first bifurcated branch of standing 
waves (SWs) from the planar base state with associated mean flow that is unforced (see 
Eq. (|56)l) and thus identically vanishes at large times. This approximation is also exact 
for the linear stability of the SWs and, in particular, for the instability threshold of this 
branch, namely the threshold amplitude (if it is finite) for the appearance of transverse 
amplitude modulations (TAMs) 5|. Furthermore, we expect that the neglected convective 
terms do not play a significant qualitative role in subsequent bifurcated branches of TAMs, 
at least near threshold. In addition this approximation is almost exact for stripes because 
convective terms can be neglected in this already explained in Sec. MI Al 

C. Two dimensional approximation 

We introduce next a (drastic) single mode approximation for the z dependence of the 
mean flow variables in Eqs. (|85 j) - (J55j) . We write horizontal velocities as 



u mo = g(z)U mo , q mo = g(z)Q mo , u ms = g(z)U ms , q mo = g(z)Q ms , (100) 
where the function g is such that 

g(z) 2 dz = l, g(0)>0, / g(z) dz > 0. (101) 

■d J -d 

This function is otherwise arbitrary and can be selected to yield the best approximation 
to the vertical velocity profiles. A reasonable choice is 



g(z) = yf2/dsm[ir{z + d)/(2d)], (102) 

which satisfies 



TV 2 



g" = -Xg, with A = — (103) 
By using this simplification Eqs. dEHJ) - dHOJ) , dH2J) and (JOHJ) reduce to 

V • U rno = 0, (104) 
d t U mo = -e 2 VQ mo + 7 (V 2 C/ mo - XU mo )/2 
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+ #l7[«(W • V) V0 + (V 2 0)V/i°) ts )' lo ], = o, (105) 
^t/ ms = -Vg ms - 7 AC/ m 72 + /3 l7 [(((V/i° • V) V0 + (V 2 0)V/i°) ts ) /ls ] z=o , (106) 
/?iQ ms - (1 - T)h ms = 0, d t h ms + p 2 V -U ms + V • {{{h°V<f>) ts ) hs ) z=0 = 0, (107) 



with 

p 1 = g(0)>0, fa 
Note that if g is given by Eq. ()102|) . then 



f g(z)dz>0. (108) 
J-d 



fa = y/2/d, p 2 = VU/tX. (109) 

The complete set of equations also includes Eq. (jB2j) and the following equations and 
boundary conditions which follow from Eqs. (|88jL (|94|) and (|9*3jl 

V 2 <P + d 2 zz <P + e 2 g"(z)(J V0) • (L/ mo + £/ ms ) = 0, (110) 

in — d < z < 0, with <9 2 = at z = —d, and 

d t h° - d z <f) + e(V • (/i°V0)) to + e 2 V • [{h mo + h ms )V<f) + p x h°{U mo + U ms )} 

+ e 2 V ■ \{h°) 2 d z V0/2] = 7 y* V 2 (d 2 0), (111) 

9,0 + + ^ 2 [(^ m ° + + (/*°) a O/2] 

+ e(| V0| 2 + |<9 2 0| 2 > to /2 + e 2 [A((7 mo + U ms ) ■ V0 + /i°«9 2 (| V0| 2 + |<9 2 0| 2 )/2] 

+ Aa(h° sin2t)*° + (1 - T)h° - TV • [Vh°/{1 + e 2 \ Vh°\ 2 ) 1/2 } + 7 <9 2 > = 0, (112) 

at z = 0. Equations ((02]) , and ()104j) - (jll2|) must be integrated with periodic boundary 
conditions in the horizontal directions, as above. Note that the mean flow does not 
contribute to the averaged energy equation at leading order (see Eq. ()35|) ). 

D. A phenomenological description 

We finally discuss the simplest possible approximation to this problem by considering 
a phenomenological model of Faraday waves that qualitatively describes its primary bi- 
furcation and secondary instabilities j^. It involves a complex order parameter ip that 
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satisfies 

d t ip = -7^ + if$ + 3i (1 + V 2 ) V/4 + (i - 7«)MV- (H3) 

A derivation of this equation is given in the Appendix. The order parameter ip is a linear 
combination of the free surface deflection and (a vertical average of) the velocity potential, 

h° = if>e- iu * + c.c., 0= -i^ + cc, or V = e~ iart (/i° + i0)/2. (114) 

Despite its simplicity, the order parameter model qualitatively describes some of the fea- 
tures of the Faraday instability: The linear dispersion relation coincides with that of the 
fluid in the limit of low viscous damping, and the model exhibits a primary bifurcation to a 
standing wave solution near threshold, which can be either sub- or supercritical at thresh- 
old depending on the wavenumber. Also, stationary solutions are in turn destabilized 
against amplitude and phase modulation instabilities. For sufficiently high supercriticali- 
ties, the solutions of Eq. ()113j) exhibit spatio temporal chaos. Consistent with the weakly 
dissipative limit we are considering in this paper, we must assume that 

7 <1, /<1, H~l. (115) 

This model has been already used by Kiyashko et al. Q] to understand how mean 
flow effects might induce rotating patterns in a Faraday wave experiment. Conjecturing 
that rotation was somehow due to the mean flow produced by surface waves, they added 
a convective term — u • to the right hand side of Eq. ()113|) . where u was a velocity 
field that was given independently of ip; thus ignoring any coupling between surface waves 
and mean flow. 

Here we shall add a similar term to the right hand side of Eq. ()113|) but, given the 
analysis above, u evolves with the surface waves according to a phenomenological equation 
with the appropriate symmetries. First we replace Eq. (jll3|) with 

d t if> = -tI) + if$ + 3i (1 + V 2 ) VV4 + (i - 7«)H 2 ^ - (3i(U mo + U ms ) • V^, (116) 

as suggested by Eqs. (jlll|) - (jll2J) . (jl!3J) and ()114|) . Note that we are not including any 
dependence on h mo and h ms , because this is beyond the scope of this phenomenological 
model. The coupling term is not conservative for general initial conditions, which is not 
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optimal. However, this term may be seen to lead to a conservative contribution at leading 
order for solutions that are linear combinations of plane waves (as in Eq. (|3H|) above). 
Since ip is intended here to only model the spatially oscillatory part of the flow, we require 
that 

(tfj) to = 0. (117) 

The contribution from the mean flow appears through XJ mo and U ms , which are the short 
and long wave components of the mean flow defined above. Their evolution is given by 
Eqs. (fTTHjl - fTTTTj) but replacing h° and by ip according to Eqs. ffTTTfl)-(fTilll 

V • U mo = 0, (118) 

3 t U mo = -e 2 VQ mo + 7 (V 2 C/ mo - XU mo )/2 + /? l7 (i(VV> • V)V^ + i(V 2 ^) Vip + c.c.) ho , 

(119) 

d t U ms = -VQ ms - 7 AC7 m 72 + /3i7(i(V^ • V)V^ + i(V 2 ^) + c.c.) hs , (120) 
(3 2 Q ms - (1 - T)h rns = 0, d t h ms + (3 2 V ■ U ms + V- (i(^V^)^ + c.c.) = 0. (121) 

We also require that (cf. Eqs. ffTH^l and (H0HJ)) 

= 16/tt 2 , A/5 2 = 4/7r. (122) 

Equations ITTT^) - ffT^TI) imply that U mo = U ms = as ft -> 0. Note also that the mean flow 
is unforced if the surface waves are standing (the phase of ip is independent of position) 
as all forcing terms in Eqs. (jll9|) - ()121j) vanish. 



IV. SECONDARY INSTABILITIES OF THE PHENOMENOLOGICAL MODEL 

In order to obtain a qualitative picture of the effects of mean flows on surface waves, 
we study in this section secondary instabilities of the base periodic solution of the order 
parameter model defined by the coupled Eqs. (jllHj) and ()118j) - (jl21j) of Sec. IIIID1 Gen- 
erally speaking, we find that mean flows couple weakly to transverse phase modulations 
and hence do not appreciably modify the zig-zag boundary. Transverse amplitude mod- 
ulations are affected by mean flows, the latter generally being destabilizing. Mean flows 
also increase the region of instability against longitudinal perturbations (Eckhaus) and, 
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more importantly, introduce a finite wave number longitudinal instability which for cer- 
tain values of the parameters can render much of the parameter space in which periodic 
solutions exist unstable. This instability branch is of oscillatory nature, and arises at 
the merging point between the branch that corresponds to long wavelength longitudinal 
modes of ip and a hydrodynamic branch which is weakly damped as k — > 0. 

The trivial solution ip = U mo = U ms = Q mo = Q ms = K ms = becomes lin- 
early unstable against a periodic perturbation of if) of wave number q for /x > /x c (g) = 



1 + [3(1 — q 2 )/^] 2 — 1, the neutral stability curve, /x is the control parameter defined 
as n = (/ — The critical mode q = 1 becomes unstable at /x = 0. 

For small /x > 0, stationary and periodic solutions exist that can be approximated by 
a single Fourier mode if> q (x) = a q cos (qx) exp (iQ q ) with 



2 _ q 2 -l- f cry 2 ± yiQf 2 (l + a 2 7 2 ) - (47 + S ai (q 2 - I)) 2 
a i~ l + (a 7 ) 2 ' 1 ' 

where the ± sign stands for sign(l — q 2 + 4a7 2 /3), and Q q satisfies sin 2Q q = (1 + 
3aa 2 /4)'y/f, cos 2Q q = |(g 2 — 1 — a 2 )/ f. Note that the bifurcation at threshold is 
subcritical if q 2 > 1 + 3«7 2 /4 ~ 1 (recall that 7 is small) and supercritical otherwise. 
This solution for the order parameter leads to vanishing driving terms in the mean flow 
Eqs. (|118JI - (|121|) : hence all mean flow variables remain zero for the base, periodic solution. 

In order to address the stability of the stationary solution (jl23j) against longitudinal 
perturbations, we introduce 

if) = yl [exp (iqx) + exp (— iqx) + a ++ exp (i(q + k)x) + a + ~ exp (i(q — k)x) 

+a~ + exp (i(fc — q)x) + a exp (—i(q + k)x)] (124) 

(where Aq = |a g exp (i0 9 )), together with the corresponding perturbations of the mean 
flow variables 

U™ s = u + exp (ikx) + c.c, (125) 
h ms = c + exp{ikx) + c.c, (126) 

and 

Q ms = d + exp (ikx) + c.c. (127) 
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where JJ mo = as seen from the incompressibility condition (jllHj) . and the definition of 
u mo = P3), which require that U™ = and (U mo ) hs = 0, respectively. 
By inserting Eq. (j!24|) into the nonlinear terms of Eqs. (jl2U|) and (jl21j) . 

=i(W • V)V^ + i(V 2 ^) + c.c, (128) 
N* =iipVip + c.c, (129) 

and retaining only the long wavelength components we obtain 

(N v x ) hs =(2id x ^J + c.c.) hs (130) 
=2\A \ 2 [{q + a ++ {q + k)e ikx + a + -{q - k)e~' lkx ){q 2 + a ++ {q + k) 2 e~ lkx 
+ a + ~(q - k) 2 e [kx ) + (-q + a- + {k - q)e ikx - a" (q + k)e-' lkx ) 
{q 2 + a- + (k - q) 2 e~ ikx + a~ (5 + A;) 2 e ifcx )] + c.c, 
(N l x ) hs =(iipd x ij + c.c.) hs (131) 
= |A | 2 [(1 + a ++ e ikx + a + -e-' lkx )(q + a ++ (q + k)e~ ikx + d + ~{q - k)e ikx ) 

+ (1 + a- + e lkx + a— e - ikx ){-q + d~ + {k - q)e~ ikx - a~ (q + k)e ikx )} + c.c, 

respectively, where N^' 1 " 1 denote the x-components of the nonlinear terms Equa- 
tions for the perturbation amplitudes u + , c + , and d + are derived by linearizing Eqs. ()12()j) . 
()121|1 with respect to all the perturbation amplitudes, after substitution of Eqs. ()13()j) 
and (j!31j) . and by extracting those terms that are proportional to e lkx : 

3 t u + = -ikd + - + 2p ll \A \ 2 [(q + £;)(2g 2 + qk){a ++ - a~) 

+ (q - k)(2q 2 - qk)(d + - - a- + )}, (132) 
p 2 d + = (1 - T)c + , (133) 
d t c + = -(3 2 iku + - iA;|A | 2 [(2g + k)(a ++ - a") + (2q - k)(a + - - a~+)]. (134) 

Finally, the governing equations for the perturbation amplitudes of the order parameter 
are obtained by linearizing Eq. (|116|) and extracting the amplitudes of the Fourier modes 
expi(±g ± k)x. For example, the governing equation for a ++ is 

d t a ++ = - 7 a++ + i/^a" + ^(1 - (q + kf)a ++ (135) 

A 4 
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+ (i - a-f)\A \ 2 (Aa ++ + 2a~~ + a + ~ + 2a~ + ) - ^\qu + . 

Similar equations result for the other three amplitudes. 

We now have a system of six first order ordinary differential equations which is linear in 
the perturbation amplitudes a , u + , and c + . The matrix of right hand side coefficients 
is denoted by A(q, k, e, . . . ), and is a function of the wavenumbers of the base solution q 
and of the perturbation k, of the control parameter e, and of the other parameters of the 
model. The base solution becomes unstable when the real part of any eigenvalue of this 
matrix becomes positive. We have numerically obtained the eigenvalues of the matrix 
A, and determined the region of stability of the base solution. Two types of instabilities 
are possible: a standard long wavelength Eckhaus instability which depends on the mean 
flow, and a finite wavenumber oscillatory instability, which is completely due to the mean 
flow. As was the case in the asymptotically exact equations for one dimensional Faraday 
waves Q], this latter instability only occurs with nonzero mean flow. Both instabilities 
will be further discussed below. 

The stability of periodic solutions against transverse amplitude and phase perturba- 
tions can be studied in a similar fashion. Given that d y ip q = in the base state with 
zero mean flow, terms involving the y components of the mean flow will be of second 
order in the amplitudes of the perturbation and hence only the components U™° and U™ s 
need to be perturbed. However, and in contrast to the case of a longitudinal perturba- 
tion, both short and long wavelength components of the mean flow need to be included. 
Furthermore, the equations for U™° and U™ s decouple at linear order, and they can be 
analyzed separately. We show next that only the U™ part modifies the transverse ampli- 
tude modulation (TAM) instability line whereas the zig-zag line is not affected by either 
component. 

We start by considering the short wavelength component of the mean flow velocity 
U™°, and introduce the following perturbation for the order parameter, 

ip = Aq [exp (iqx) + exp (— iqx) + a ++ exp (i(qx + ky)) + a + ~ exp (i(qx — ky))+ 

+a~ + exp (i(— qx + ky)) + a exp (— i(qx + ky))\ , (136) 

U™° = v ++ exp {i(2qx + ky)) + v + ~ exp (i(2qx - ky)) + c.c, (137) 
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and 

quo = p ++ exp ^ 2qx + k y ))+ p+- eX p (i(2gx - ky)) + c.c. (138) 

The amplitude of the mode e l2qx in the x-component of the nonlinear forcing term of Eq. 
(HH is 

(N v x ) h ° = ([d x ^ + d y ^d xy ^+ (d B J+dffid^\ + c.c.} h0 (139) 
= \A \ 2 qk 2 [(a- + - a + -)e- [ky + (a~ - a ++ )e iky }e i2qx + c.c., (140) 

where in Eq. (I14UJ1 only terms linear in a ±± and a ±=l= have been kept. The pressure Q mo 
is calculated by taking the divergence of Eq. (|119|) . 

= -V 2 Q mo + /3^(d x N: + d y N v y ), (141) 

and is thus eliminated by using 

= (Aq 2 + k 2 )p ++ + A7i2g|A | 2 gfc 2 (a~ - a ++ ), (142) 

and a similar equation for p + ~ (here we have used the fact that one obtains (Ny) ho = at 
linear order). We then derive the following linear system of equations for the perturbation 
amplitudes, 

d t a ++ = - 7 a ++ + i/^a— + -(1 - q 2 - k 2 )a ++ 

Aq 4 

+ (i - a7)|A,| 2 (4a ++ + 2a~ + a + ~ + 2a~ + ) + P^qv^. (143) 

d t v ++ = - l{Aq 2 + k 2 + X)v ++ + p a \A \ 2 -^—(a~ - a++), (144) 
2 4q z + At 

with similar equations for a + ~, a~ + , a ++ , and v + ~ . 

A transverse amplitude modulation (TAM) is defined by the linear combinations b\ = 

a ++ + a + ~ + a' + + a~ and v 1 = ^s(v ++ + v + ~). From Eqs. and (jT4H) we find a 

closed system 

A Q * Q 

OA = -7&1 + U^h + y(l - q 2 - k 2 )h + -« 2 (i - a 7 )(26i + h) - 2(3 iqVl , (145) 
/±0 4 4 

^ = -|(4g 2 + fc 2 + A),, - ^^f^fe, (146) 

The relevance of mean flows on this perturbation is demonstrated in Fig. showing the 
stability boundary in the plane (/3i,/i) at fixed wavenumber q and model parameters. 
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With increasing mean flow coupling (3i, the region of stability of the base solution (region 
above the dashed line in the figure) decreases. 

Alternatively, a transverse phase modulation (zig-zag) is given by b 2 = a ++ — a + ~ — 
a~ + + a and v 2 = Q(v ++ — v + ~). We find in this case, 

A Q * 1 

d t b 2 = - 7 &2 + if-^h + -(1 - q 2 - k 2 )b 2 + -a 2 q (i - aj)(2b 2 + b 2 ) - 2^qv 2 , (147) 

-Tig 4 4 

d t v 2 = + k 2 + X) V2 - ^-^-^ (148) 

From the eigenvalue equation of this system, it can be shown that modifications to the 
eigenvalues due to v 2 appear at higher order in k. Since the zig-zag instability occurs in 
the limit fc^Owe find the effect of mean flow to be negligible in this case. 

A similar analysis has been performed for the long wavelength component of the mean 
flow with perturbations of the form U™ s = u + exp (iky) + c.c, h ms = c + exp (iky) + c.c, 
and Q ms = d + exp (iky) + c.c. The order parameter is again given by Eq. ()136|) . As 
shown in Fig. ^ the contribution of the long wavelength component of the mean flow to 
transverse modulations is negligible. 

We now turn to a summary of our numerical results about the base periodic solution. 
Figure El shows the various stability boundaries for the special cases of (3\ = (no mean 
flow) and fli = 0.5. Other values of the parameters used are 7 = 0.1, a = 0.5, and 
T = 0.8. The values of /3 2 and A depend on /3i according to Eq. ()122j) . Except for 
a, these parameter values correspond approximately to the values for the low-viscosity 

n 

experiments of Kudrolli and Gollub described in [5]. For instance, typical experimental 
values of (3\ = (where d is the dimensionless height of the layer) are between 0.5 

and 1.2. Figure El includes the neutral stability curve of the basic periodic solution and, 
since the primary bifurcation is subcritical for q > 1, we have included the saddle node 
curve where the periodic solution bifurcates (a subcritical bifurcation for q > 1 as 7 — ► 
has also been found in a direct numerical solution of the governing fluid equations in two 
dimensions The case /3i = is shown as a reference, and it agrees with the results 

of H- 

The range of base solutions that is stable against all perturbations considered here 
(Eckhaus, TAM, and zig-zag) is a small region close to threshold at fi = between the 
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TAM and zig-zag lines. Periodic solutions are stable against transverse perturbations 
below the dashed-dotted line in the figure (zig-zag, denoted "Z"), and above the dashed 
line (TAM). Eckhaus perturbations have negative growth rate below the dotted line. 
We observe that with increasing /3i both Eckhaus and TAM curves are shifted so that 
larger regions in the (/i, q) space become destabilized with respect to TAM or Eckhaus 
perturbations. As discussed above, the zig-zag line is not affected by the mean flow. 

We finally discuss a new oscillatory instability against longitudinal perturbations which 
is absent for /3i = 0. The oscillatory nature of the instability is demonstrated in Fig. |3] 
that shows the real and imaginary parts of the corresponding critical eigenvalue branch 
as a function of the wavenumber of the perturbation k for fixed q, fi, and other model 
parameters. The imaginary part of the eigenvalue a is not zero at the point in which 
3?(cr) = 0. This figure also shows that the instability occurs at small but finite wavenumber 
k, fact that has been confirmed by calculating the wave number with largest growth rate 
both slightly above and below the instability threshold at fi = 0.1155. 

The inset in Fig. El illustrates the origin of the oscillatory instability. In the limit of 
k — > there are two distinct eigenvalue branches that have a small and negative real 
part. The upper branch is marginal at k = and is related to the translational symmetry 
broken by the base state if) q {x). On the other hand, the mean flow velocity vanishes in the 
base state, originating the lower branch which is weakly damped at k = 0. The damping 
rate of the relevant mode for longitudinal perturbations (U™ s ) is 7A/2 as can be seen 
from Eq. ()120|) . As k increases, the two (real) eigenvalue branches merge leading to a 
complex conjugated pair and to an oscillatory instability. 

If the eigenvalue problem for the perturbations of Eckhaus type is expressed by the real 
and imaginary parts of the coefficients a ±± , , u + , and c + two pairs of complex conjugated 
eigenvalues cross the imaginary axis at the instability point. Therefore, the effective di- 
mensionality of the critical subspace is four. This can be understood from the symmetry 
group of the system of Eqs. ()132J) - (jl36|) . These equations are invariant under a spatial 
reflection (exchanging a ++ with a , a + ~ with a h , and u + with — u + ) and rotation (act- 
ing on a ±+ and u + as multiplication by exp (iO) and on a ±_ by exp (— iB), where O is 
arbitrary) that derive from the symmetry of the original equations. The corresponding 
symmetry group 0(2) has two-dimensional irreducible representations which in turn re- 
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quires each eigenvalue to occur twice. At threshold, the eigenspace of the linear system 
is spanned by four linearly independent eigenvectors, two for each pair of the complex 
conjugate eigenvalues. Symmetric bifurcation theory shows that for generic bifurcations 
with 0(2) symmetry the nonlinear solution branches are standing and traveling waves 
j^l . Examples of the temporal evolution of eigenvectors of both types are plotted in Fig. 
@] In practice, the specific values of the parameters determine the way the standing and 
traveling waves bifurcate If both are supercritical, the one with larger amplitude 

would be stable, the other one unstable. 

We finally show our results concerning the location of the oscillatory instability bound- 
ary as a function of q for two values of fi\ in Fig. Stable regions are located to the right 
of the plotted curves. As expected, the stability boundary moves toward the Eckhaus 
line with decreasing /3i, presumably merging with it for f3\ — > 0. Note that for f3\ = 0.5 
the unstable region covers most of the region of existence of the base states except for a 
narrow stripe close to the saddle node bifurcation. 



V. DISCUSSION 



The coupled system of equations describing fast surface oscillations and slowly evolving 
mean flows in three dimensions has been derived. Mean flows are forced by the surface 
waves by various mechanisms. For a particular choice of geometry and limit of parameters 
which are relevant to recent experiments, we have shown that two contributions appear 
at the appropriate order in the multiple scale expansion: A viscous streaming flow forced 
at the free surface with components of both similar and larger scale compared to the 
length scale of the surface waves, and a long wavelength component originating from slow 
distortions of the surface elevation that exists even in the absence of viscosity. 

The analysis presented has illustrated the importance of mean flows in the Faraday 
wave system for small viscous damping by determining the stability boundaries of the base 
pattern of standing waves against long wavelength perturbations. Since the full system 
of surface wave/mean flow equations is quite involved, we have instead carried out the 
stability analysis of a phenomenological order parameter equation, similar in spirit to the 
Swift-Hohenberg model of Rayleigh-Benard convection. In addition, we have limited the 
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analysis to the simplest regular pattern consisting of stripes. Mean flows are induced by 
perturbation of the stripes, and their coupling to the order parameter equation affects 
the stability of stripe solutions. We have found that the mean flows generally destabilize 
the base solution. The strongest coupling, and hence the strongest destabilization, occurs 
for longitudinal or Eckhaus perturbations. Furthermore, mean flows introduce a new 
oscillatory instability which for small nonlinear damping in the phenomenological model 
renders all stripe patterns unstable. A weaker effect has been found for finite wavelength 
transverse amplitude modulations which largely couple only to the short wave part of the 
mean flow. 

Our first remark concerning experiments follows from the existence of a longitudinal 
oscillatory instability. Within the phenomenological model, all stripe solutions which are 
stable in the absence of mean flows are unstable against longitudinal oscillatory perturba- 
tions for sufficiently large coupling parameter f3\. One would then expect time dependent 
behavior at onset. The eigenvectors corresponding to the oscillatory instability (Fig. QJ 
show that the associated mean flow consists of large scale rolls with their axis oriented 
parallel to the surface. At the surface it advects the waves leading to compression and di- 
lation of waves similar to Eckhaus perturbations, but in the form of traveling or standing 
waves. A numerical solution of the coupled order parameter mean flow equations shows 
that the compression not only leads to a decrease in wave amplitude, but can also result 
in a complex cycle including the annihilation of stripes, possibly due to a different in- 
stability triggered by the compression. We also anticipate novel phenomena arising from 
mean flows if one considers the slow dynamics of defects in the wave pattern. Defects 
as local perturbations of a regular pattern drive mean flows which in turn affect defect 
motion. 
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APPENDIX A: ORDER PARAMETER MODEL 



We briefly summarize in this Appendix the derivation of the phenomenological model 
given by Eq. ()113|) of Sec. IIIIDl first introduced in ref. [3. We follow the description 
originally introduced by Zakharov [3^, and later of Crawford, Saffman and Yuen 3^| 
in their study of the nonlinear evolution of deep water gravity waves in an inviscid, 
incompressible, and irrotational fluid. Their analysis can be straightforwardly extended to 
a parametrically driven fluid, and linear viscous damping is added in a phenomenological 
way. All the variables used in this Appendix are assumed to be dimensional quantities. 

The governing equation for the inviscid fluid is, 

(V 2 + <9 2 )0 = O, -00 < z < /i(x,i), (Al) 
with boundary conditions at the free surface z = /i(x, t), 

d t h + Vh-V(f) = d z (f), (A2) 

M + \ (V0) 2 + \ {d^f + (g + g z (t)) h = -V- ] , (A3) 

with a the interfacial tension and p the density of the fluid that is being vibrated with 
acceleration (— g — g z {t)) in the z direction. It is well known that this problem admits a 
Hamiltonian formulation with Hamiltonian 

H= l - Jdxp 1 * dz (V0) 2 + [d z <pf + l(g + g z (t)) h 2 + - (V 1 + (Vhf - 1^ 

P ' (A4) 

where the velocity potential further satisfies the boundary condition d z (f) = as z — > —00. 
The canonically conjugate variables are the surface deflection h(pc, t) and the velocity 
potential on the surface s (x, t) = </>(x, z = h(pc)). Phenomenological damping can be in- 
troduced by considering a dissipation function Q(h(x, t), s (x, £)). The resulting canonical 
equations of motion are, 

a < hM =wh' (A5) 

<9 t s (x, t) = "^^y + Q(Hx, t), s (x, t)). (A6) 
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The functional Q determines the rate of viscous dissipation in Eqs. (|A5|) - (|A6|) so that 

J dxQ(h(x,t),<f) s (x,t))d t h(x,t). (A7) 



dH dH 
~dt dt 



Of course, Q = corresponds to the inviscid limit. 

The case of a fluid of low viscosity has been treated by assuming that energy dissipation 
is dominated by potential flow in the bulk [29] . The functional Q can then be determined 
by equating the rate of dissipation in Eq. (jA7|) to the rate of energy dissipation due to 
potential flow, 

J dxQ(/i(x,£),0 s (x,t))<9 t /i(x,t) =-v J dxj ( ' } ckV 2 (V0) 2 . (A8) 

This equation has been used to determine Q order by order in an expansion in the surface 
wave steepness 0, Q| • To the order relevant here, one finds, 

Q(k, t) = -4vk 2 (f) s (k, t) + nonlinear terms, (A9) 



where Q(k, t) is the Fourier transform of Q. As discussed in refs. |10|, [ll|, this approxi- 
mation yields, by construction, the correct rate of energy dissipation at linear order, but 
not the correct equations of motion even at this order. In particular, it overestimates by 
a factor of two the damping force in Eq. (|A6J) . and omits wave rectification in Eq. (JA5j) 
that arises from the rotational component of the flow in a thin boundary layer adjacent 
to the free surface. 

Following Zakharov j^], we define a complex field 



KM) = /fU( M ) +i ^-^( M ), (AlO) 
where h(k, t) and <j) s (k, t) are the two dimensional Fourier transforms of /i(x, t) and <p s (x, t) 



respectively, and u(k) = g k + ak 3 /p is the inviscid dispersion relation. In terms of 
this new variable, the Hamiltonian system ()A5j) - (jA6j) can be written as, 



9 ' 6(k - t) = - ! .TO + ! V^ (k ' i) - (AU) 

Equation (jAll|) is now expanded in a power series of b. We confine ourselves here to linear 
terms in b, as nonlinear terms will be added phenomenologically. However, explicit forms 
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of cubic terms in b have been obtained |42| both for the present case of an expansion 
around the inviscid solution, and also for the Linear Damping Quasi Potential Equations 



of ref. 



y 



By expanding Eq. (jAll|) in power series of b, we find, 

d t b(k, t)+2uk 2 [b(k, t) - 6(-k, t)] +iu(k)b(k, t)+^^r *) + 6(-k, t)] +JVX(6(k, t)) = 0, 

(A12) 

where Af£(b(k, £)) stands for terms nonlinear in the amplitudes b. 

If the driving acceleration is given by g z (t) = a cos fit, only amplitudes with wavenum- 
ber close to the critical wavenumber ko are excited near onset, with frequency close to the 
resonant frequency u(ko) = 0/2. We introduce a conventional multiple scale expansion 
near onset, but choose to do so so in a manner that will preserve the rotational invari- 
ance of the original governing equations. We further assume the following scalings for the 
damping and driving terms: 7' = 2uk 2 = e 2 7o and /' = k a/ (Auj(ko)) = e 2 fo, where e is a 
small expansion parameter, and both 70 and fo are 0(1) quantities. We also expand 

6(k, t) = eB(k, T u T 2 )e-^ fc °)' + e 2 6 2 (k, t) + e 3 6 3 (k, t) + . . . (A13) 

with T\ = et and T 2 = e 2 t. The slow time scale T\ corresponds to the time scale of 
translation of a wave packet, whereas T 2 is the scale of change in the modulation of 
the wave packet. These two time scales are consistent with an expansion of the inviscid 
dispersion relation u(k) = uj(ko) + euo' + e 2 u" + . . . for modes near the critical wavenumber 
ko- Substitution of Eq. (|A13|) into Eq. (|A12|) shows that the equation is identically 
satisfied at 0(e). At 0(e 2 ) we obtain the following solvability condition, 

f) Ft 

^— = -iuj'B. (A14) 

The solvability condition at order 0(e 3 ) is 

— = -7oS(k) -if B(-k) -iu"B(k) +NC[B}. (A15) 
01 2 

with a known nonlinear functional MC[B]. We now combine the two solvability conditions 
by writing A(k) = eB(k) and d t A = e 2 dx 1 B + e 3 dT 2 B, and find, 
8 A 

^— = -7 , A(k, t) - if'A(-k, t) - % [eJ + e 2 oo") A(k, t) + MC[A] (A16) 
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Hence the slow evolution near onset given by Eq. (jA16|) is the same as that of the original 
set of inviscid equations supplemented by phenomenological linear damping. 
From the inviscid dispersion relation, we find, 

ecu' + eV = c x (k 2 - kl) + c 2 (k 2 - klf + 0{{k - A; ) 3 ), (A17) 

with 

9o 3ak 

Cl = IhMkoj + ' (A18) 

and 

1 / 3afe {g + 3ak 2 /p) 2 g \ 
° 2 4k 2 u(ko) V 4p 8u(k ) 2 4k J- { ] 

We derive next a real space order parameter model from Eq. (|A16|) . Three simplifica- 
tions are necessary. First, the nonlinear functional in Eq. (|A16|) does not have a closed 
form representation in real space. As has been done in other systems (cf. Rayleigh-Benard 
convection, j^J), we introduce phenomenological functional forms for this term. In doing 
so we artificially determine the symmetry of the bifurcating pattern at onset, but more 
importantly in the case of Faraday waves , w e sidestep the issue of the origin of nonlinear 
damping and saturation of the waves ^H^jJ- I n ^ ne simplest possible case, the nonlinear 
term in Eq. (|A16J) is approximated by an imaginary constant iR'. Second, it is also known 
that linear damping is not sufficient to produce wave saturation in this system 40]. We 
introduce a phenomenological nonlinear damping coefficient 07', where a is a constant 
assumed to be of order 1. 

We finally define a complex order parameter field tp(pc,t) as the inverse Fourier trans- 
form of A(k), and find from Eq. (|A16|) . 

d t tp = -7V + if i) + id + V 2 ) ij - ic 2 (V 2 + kl) 2 if> + (-aV + iiO||V||V- ( A2 °) 

We now choose uj(ko) = ft/2 as the unit of time, l/fco as the unit of length, and further 
define / = f'/uj(k ) = k a/4:Uj(ko) 2 , 7 = ^ /uo(k ) = 2vkl/uj(k ), and R = R'/u(k ). By 
choosing 1 / \^R as the order parameter scale, Eq. (jll3|) results in the capillary wave limit. 
The positive sign of the imaginary part of the nonlinear coefficient (— aj+i) in Eq. (jll3|) is 
chosen to represent capillary waves. In the opposite limit of gravity waves, the imaginary 
part of this coefficient has to be negative. Note that as a third simplification we have 
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eliminated the term i (l + V 2 ) 2 ip in Eq. (J113j) as this term together with i (l + V 2 ) ip 
leads to two different wavenumbers becoming critical at threshold, and unwanted feature 
for us. 

The effect of this third simplification can be further understood by comparing the 
amplitude equation on the model and that of the inviscid fluid. By introducing a multiple 
scale expansion of the form 

±N 

?P = 5J2 %( X > *> T i> T 2 )^ + 5 2 ^ 2 + 6 3 tfj 3 . . . (A21) 
j=±i 

with 5 a small bookkeeping parameter, and X = £x, 7\ = 5t, and T2 = S 2 t we find that 
up to order 0{5 3 ) 

dtdj = -7o 3 - + ifa-j - ^ (kj ■ V) ^ + j V 2 aj + 

+ (— «7 + z) j |%| 2 aj + 2 |a/| 2 aj + 2 a;a_/a_j J (A22) 

by following the same expansion procedure outlined above. The terms linear in the ampli- 
tudes are the same as the corresponding terms in the amplitude equation derived directly 
from the inviscid equations except for an additional term (kj • V) 2 aj which is missing in 
Eq. (|A22j) . This is a direct consequence of having eliminated the term + V 2 ) 2, in 
the phenomenological model. 
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FIG. 1: This figure illustrates the weak destabilization of the base solution against Transverse 
Amplitude Modulations due to the coupling to mean flows. Shown are the critical values for 
instability [i C r%t arising from either the long wavelength component of the mean flow velocity 
alone ("ms") or the short wavelength component alone ("mo") as a function of the coupling 
parameter /3i. The base solution is stable above the corresponding lines. The figure shows that 
the "ms" component of the flow does not appreciably modify the stability threshold, whereas 
the effect of the "mo" component is to weakly destabilize the base solution with increasing (3\. 
The wavenumber of the base solution is q = 1.04, and we have used 7 = 0.1, a = 0.5, and 
r = 0.8. The parameter (3\ spans the range between no mean flow and the approximate value 
that corresponds to the weak viscosity but shallow layer experiments of ref. 
/?2 and A depend on (3\ according to Eq. (|122|) . 
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FIG. 2: Stability diagram for (a) the order parameter model of Eq. without mean flow, 

and (b) with f}\ = 0.5. Other values of the parameters used are 7 = 0.1, a = 0.5, and V = 0.8. 
We show the Eckhaus line (dotted) , TAM line (dashed) , and the zig-zag line (dash-dotted) . We 
also show the neutral stability curve of the primary instability (solid line denoted by N), and a 
saddle node bifurcation (thick solid line marked S). Only the left branch of the Eckhaus line is 
shown emanating from {q = = 0). This line is parabolic near the critical point, but quickly 
bends to the right as shown in the figure. Hence the region of stability of the base solution 
against an Eckhaus instability is the region below the dotted line. Comparison of (a) and (b) 
shows that the mean flow decreases the regions of stability against Eckhaus modulations and to 
a small extent against transverse amplitude modulations. In (a) the region of stability against 
all perturbations is shown by the gray area. This region is not indicated in (b) since these 
solutions are unstable with respect to the oscillatory instability. 

FIG. 3: Real (left) and imaginary (right) part of the largest eigenvalue for f3\ = 0.2, \x = 0.1318 
and q = 1.07. The inset shows the region near k = 0. 

FIG. 4: Critical eigenvectors corresponding to the oscillatory instability at f5\ = 0.2, q = 1.0, 
/i = 0.1318, and k = 0.04. We show the temporal evolution over its own period (vertical 
axis) to illustrate the nature of the critical modes. The eigenvectors are represented by the 
x component of u ms (a,c) and the amplitude of ip (b,d). The evolution is obtained from the 
two critical eigenvectors Vi of the complex matrix A(q, k, . . . ) with eigenvalues Left and 

right traveling (a,b) waves are given by the evolution of v\e lu)t and V2e~ lujt , respectively. A 
standing wave is obtained (c,d) by superposition of the evolution from both eigenvectors scaled 
to equal amplitude. The actual solution for tp is obtained from the superposition of the plotted 
eigenvectors with the base solution, see Eq. ()124l) . 

FIG. 5: Stability boundaries of the oscillatory instability for three values of (3\\ 0.05 (dotted 
line), 0.2 (dashed line), and 0.5 (solid line). The thick solid line indicates the saddle node. 
Periodic solutions are unstable to the left of the curves. Comparison with Fig. |2] shows, that all 
periodic solutions are unstable and the pattern is expected to be time-dependent. 
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